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Abstract 

Continuing previous work on the 3PN-accurate gravitational wave generation from point particle 

■ binaries, we obtain the binary's 3PN mass-type quadrupole and dipole moments for general (not 
necessarily circular) orbits in harmonic coordinates. The final expressions are given in terms of 
their "core" parts, resulting from the application of the pure Hadamard-Schwartz (pHS) self-field 
regularization scheme, and augmented by an "ambiguous" part. In the case of the 3PN quadrupole 
we find three ambiguity parameters, £, k and £, but only one for the 3PN dipole, in the form of 
the particular combination £ + k. Requiring that the dipole moment agree with the center-of-mass 

■ position deduced from the 3PN equations of motion in harmonic coordinates yields the relation 
0^ . £ + k = —9871/9240. Our results will form the basis of the complete calculation of the 3PN 

\ radiation field of compact binaries by means of dimensional regularization. 
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I. INTRODUCTION 



The present paper is the continuation of previous work [1] 1 on the generation of grav- 
itational waves by inspiralling compact binaries, viz. binary systems of neutron stars or 
black holes whose orbit adiabatically spirals in by emission of gravitational radiation. The 
adiabatic inspiral takes place right before the final plunge and merger of the two compact 
objects, to (presumably) form a single black hole which will settle down, after emission of 
its quasi-normal modes, into a stationary configuration. Inspiralling compact binaries will 
almost surely be detected by large scale laser interferometric gravitational wave observato- 
ries like VIRGO and LIGO. Recent estimates of the rate of coalescences of neutron stars are 
very promising 0. 

It is by now well established (see e.g. Refs. 0, IE El 0, GU I3|) that, in order to 



predict in a useful way the gravitational radiation emitted by inspiralling compact binaries, 
general relativity must be developed to high post-Newtonian (PN) order, probably up to 
the 3PN or even the 3.5PN level. 2 Correlatively it has been realized that the crucial point 
for detecting and deciphering the gravitational waves, is to accurately take into account the 
gravitational interaction between the compact objects — responsible for the binary's orbital 
dynamics and wave emission. Arguments within general relativity [12| show that a good 
modelisation is by point particles, characterized only by mass parameters mi and m-i (we 
neglect the intrinsic spins of the compact objects). It makes sense to implement a model 
of point particles within the PN approximation, provided that a process of regularization is 
used for dealing with the infinite self-field of the point particles. The regularization should 
hopefully be followed by a renormalization. 



In paper I we adopted as self-field regularization the Hadamard regularization [13J, |1J, [11 
augmented by a prescription for adding a few arbitrary unknown "ambiguity parameters" , 
accounting for the incompleteness of the Hadamard regularization when evaluating certain 
divergent integrals occuring at the 3PN order. We found that the 3PN mass quadrupole 
moment of point-particle binaries is complete up to three ambiguity parameters, denoted £, 
k and (, which would typically be some rational fractions and could take, within Hadamard's 
regularization, any numerical values. [The quadrupole moment is the only one to be com- 
puted with full 3PN accuracy, thus it contains most of the difficult non-linear integrals, and 
all the ambiguities associated with the Hadamard regularization.] The gravitational-wave 
flux, which is a crucial quantity to be predicted because it drives the binary's orbital phase 
evolution, has then been found to be complete, in the case of circular orbits, up to a single 
combination of the latter ambiguity parameters, given by 9 = £ + 2k + (. Of course the 
ambiguity parameters do not affect the test mass limit of the result of paper I, which is 
found to be in perfect agreement with the result of linear black-hole perturbations in this 



round to be m per : 
limit HE Bill 



The parameters £, k and ( represent the analogues, for the case of the gravitational 
wave field (more precisely the mass quadrupole moment), of similar parameters which were 
originally introduced in the problem of the equations of motion of point particle binaries 
at the 3PN order More precisely, £ and k are the coefficients of some 



static terms, independent of the particle's velocities but depending on their accelerations, 



1 Henceforth the Ref. [l| will be referred to as Paper I. 

2 As usual the roPN order refers to the terms of order 1/c 2 " in the waveform or energy flux, relatively to 
the lowest-order (Newtonian) approximation described by the Einstein quadrupole formula. 



2 



which can be viewed as some analogues of the "static" ambiguity constant u s in the 3PN 
ADM-Hamiltonian (lil 2oT ]. which is itself equivalent to the parameter A entering the 3PN 
equations of motion in harmonic coordinates 0, 2^] . On the other hand, ( is the coef- 
ficient of a particular velocity-dependent term in the mass quadrupole, and is similar to 
the "kinetic" ambiguity uj^ in the Hamiltonian 0, 2(|, which has no counterpart in the 



equations of motion since the velocity terms were unambiguously determined there (21. 22l . 
The work (2ll . 22l | used an improved version of the Hadamard regularization called the 
extended-Hadamard regularization, based on a theory of pseudo-functions and generalized 



tneory 

US 



distributional derivatives, and defined in Refs. 

The ambiguity parameters in the binary's local dynamics (Hamiltonian and/or equations 
of motion) have been resolved. For the kinetic ambiguity we have uik = 41/24, which 
follows from the requirement of invariance under global Poincare transformations [2jJ, |2j| . 
On the other hand the static ambiguity has been fixed using a powerful argument from 
dimensional regularization, i.e. computing the binary's dynamics in d = 3 + e spatial 
dimensions and considering the limit where e — > 0, which led to to s = [26] or, equivalently, 
to A = —1987/3080 27]. The same result has also been achieved in Refs. |28l l29l. I30I l3ll| by 
means of a surface integral approach to the equations of motion of compact objects (i.e., a 
la Einstein- Infeld- Hoffmann), successfully implemented at the 3PN order. 

Summarizing, the 3PN equations of motion have been completed in essentially two steps: 
the first one consists of using Hadamard's regularization and permits the computation of 
most of the terms but a few; the second step is to apply dimensional regularization in order 
to fix the value of the few parameters left undetermined. In the present situation it is not 
possible to compute the equations of motion in the general <i-dimensional case, but only in 



the limit where e = d — 3 — > |26l . l27l |. In Refs. |26l . |27| one computes the difference be- 
tween the dimensional and Hadamard regularizations, and it is this "difference" , specifically 
due to the existence of poles in d dimensions (proportional to 1/e), which corresponds to 
the ambiguities in Hadamard's regularization. Actually the latter difference has to be de- 
fined with respect to a particular Hadamard-type regularization of integrals called the "pure 
Hadamard- Schwartz" (pHS) regularization, following the terminology and definition of Ref. 
(2?}. The pHS regularization consists of the standard notion of Hadamard's partie finie of di- 
vergent integrals, together with a minimal treatment of the compact-support (or "contact" ) 
terms, and the use of Schwartz distributional derivatives [Til]. The result of dimensional reg- 
ularization is given by the sum of the pHS regularization and of the "difference" containing 
the poles proportional to \je. 

In the present approach, for the 3PN wave generation we basically follow the same two- 
step strategy as for the 3PN equations of motion, namely: 

(i) To obtain the expression of the mass quadrupole moment at 3PN order, as regularized 
by means of the pHS regularization; 

(ii) To add to the pHS result the difference between the dimensional regularization and 
the pHS one, which as we said above is due to the presence of poles at the 3PN order. 

Imposing then that the result of dimensional regularization is equivalent to the result of 
(the pHS variant of) Hadamard's regularization and augmented by appropriate ambiguity 
parameters, will then uniquely determine the ambiguity parameters. A summary of the 
above calculations leading to the following unique values for the ambiguity parameters £, k 
and (: 

9871 7 . , 

* = _ 9240' K = °' C = ~33' (L1) 
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has been provided in Ref . |32| . The technical details of the above calculations are now given 
in a series of three papers of which the present one is the first. Indeed, the present paper is 
devoted to the calculation of the pHS regularization of the 3PN quadrupole moment, item (i) 
above. The computation of the part associated with the difference between the dimensional 
and pHS regularizations, cf. item (ii), will be provided in the third paper of this series j33|. 
In the second paper |34| the value of ( has been confirmed by a different approach, based 
on the multipole moments of a boosted Schwarzschild solution. 

It should be emphasized that the values (jl.lj) represent the end result of dimensional 
regularization, obtained by means of the sum of the steps (i) and (ii). However, we shall be 
able to obtain below [see Eq. (jl.2|) and Section IV Bj an independent confirmation, within 
Hadamard's regularization, of the value of the particular combination of parameters £ + k. 
Moreover, the fact that k = has been checke d by a diagrammatic reasoning in Ref. j33|. 
Since as we said ( has also been computed in [34| by a different method, we see that the 
present paper and the works |33|, |34| altogether provide a check, independent of dimensional 
regularization, for all the parameters (jl.lj) . 

In our previous work (paper I), the 3PN mass quadrupole moment was regularized by 
means of some "hybrid" Hadamard-type regularization, instead of the pHS one, and the 
ambiguity parameters £, k and ( were defined with respect to that regularization. In the 
present paper, since we shall perform a different computation, based on the specific pHS 
regularization, we shall have to introduce some new ambiguity parameters. Since we do not 
want to change the definition of £, k and (, we shall perform some numerical shifts of the 
values of £, k and (, in order to take into account the different reference points for their 
definition: hybrid regularization in paper I, vs. pHS regularization in the present paper. 

The present investigation will also extend and improve the analysis of paper I in two 
important ways. First we shall use a better formulation of the multipole moments at the 
3PN order, in terms of a set of retarded elementary potentials, instead of the "instantaneous" 
versions of these potentials as was done in paper I. The retarded potentials are the same as 
in our computation of the equations of motion in harmonic coordinates j^; their use will 
appreciably simplify the present work. Secondly we shall generalize paper I to the case of 
arbitrary orbits, not necessarily circular. Circular orbits are in principle sufficient to describe 
the inspiralling compact binaries, but the general non-circular case will be mandatory when 
we want to obtain the values of £, n and ( separately, and it is also important for a check 
of the over-all consistency of our calculation. 

Besides the 3PN mass quadrupole moment of point particle binaries we compute also their 
3PN mass dipole moment. The mass dipole is interesting because it is a conserved quantity 
(or it varies linearly with time), which is already known from the conservative part of the 
binary's local 3PN equations of motion in harmonic coordinates. Namely the dipole moment 
is nothing but the integral of the center- of-mass position associated with the invariance of 
the equations of motion under the Poincare group. It has been computed from the binary's 
Lagrangian in harmonic coordinates at 3PN order in Ref. 3^, 36| . In fact what we shall do 



in the present paper is to impose the equivalence between the 3PN dipole moment and the 
3PN center-of-mass vector position, and we shall prove that this requirement fixes uniquely 
one, but only one, combination of the ambiguity parameters, viz. 

9871 , . 

f + K = -9240' (L2) 

working solely within Hadamard regularization. This result is perfectly consistent with the 
complete result provided by dimensional regularization [32[ and recalled in Eq. (jl.lj) . We 
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view this agreement as an important check of the correctness of both the present calculation 
and the one of Ref . j32] • 

The paper is organized as follows. In Section|n]we review our definitions of the multipole 
moments of an isolated extended source in the PN approximation, both time-vaying moments 
(having i > 2) and static ones (£ = 0, 1). In Section ITTrj we give the explicit expression of 
the mass-type moments in terms of a set of retarded elementary potentials up to 3PN order. 
The pure-Hadamard-Schwartz (pHS) regularization scheme is then reviewed in Section I1V1 
where we comment on the various types of terms encountered in the calculation, and we 
detail our practical way to perform the partie finie of three dimensional non-compact support 
spatial integrals. Our final results for both the 3PN quadrupole and dipole moments, and 
our derivation of Eq. (jl.2|) . are presented in Section|3 The formula for the 3PN quadrupole 
moment in a general frame turned out to be too long to be published, therefore we choose to 
present it in the frame of the center of mass (but for general orbits): see Eqs. (|5.9|) - (j5.1U|) 
below. 



II. MULTIPOLE DECOMPOSITION OF THE EXTERIOR FIELD 

In this Section we provide an account of the relevant notion of multipole moments of 
a general isolated gravitational wave source. By definition the moments parametrize the 
linearized approximation in a post-Minkowskian expansion scheme for the gravitational field 
in the external (vacuum) domain of the source. Their explicit expressions in terms of the 
source's physical parameters (matter stress-energy tensor T^ u ) have been found in the case of 
wst-Newtonian source by using a variant of the theory of matched asymptotic expansions 
3^| . The matching relates the exterior field of the source, as obtained from a 



multipolar-post-Minkowskian expansion of the external field |4£|, to the inner field of the 
post-Newtonian extended source, as iterated in the standard PN way. 



A. External solution of the field equations 

The Einstein field equations are cast into "relaxed" form by means of the condition of 
harmonic (or De Donder) coordinates. Denoting the fundamental gravitational field variable 
by W = \J—g g^ v — iff" , 3 this means that 

d v h" v = 0. (2.1) 

Under the harmonic coordinate conditions the field equations take the form of non-linear 
wave equations 

UhT = t^L T ^ j (2.2) 

in which □ = f] pcr d pcr denotes the standard flat space-time d'Alembertian operator. The 
right-hand-side of Eq. (|2.2|) is made of the total (matter plus gravitation) pseudo stress- 
energy tensor in harmonic coordinates given by 

r 4 

T f" = \ g \T^ + -A"" \h, dh, d 2 h] , (2.3) 

16nG L J 



3 Here, g^ v denotes the inverse of the usual covariant metric g M „; g is the determinant of g^ v : g — det(g fJiU ); 
and rf ,v is an auxiliary Minkowskian metric in Minkowskian coordinates: ff" v = diag(— 1, 1, 1, 1). 



5 



where T^ v is the stress-energy tensor of the matter fields, and A pv represents the gravitational 
source term which is given by a complicated non-linear, quadratic at least, functional of h pa 
and its first and second space-time derivatives. Equation (2.3) in paper I gives the explicit 
expression of A^ v . In the following we shall assume that the support of T M1/ is spatially 
compact. In our formalism the conservation of the pseudo tensor, 

d v r^ = , (2.4) 

is the consequence of the harmonic-coordinate condition (j2.1|) . 

Let the calligraphic letter M. denote the operation of taking the multipole expansion, 
so that M.(h^ u ) represents the multipole expansion of the external gravitational field - 
a solution of the vacuum field equations valid outside the compact support of the matter 
tensor T^ u . Similarly denotes the multipole expansion of the gravitational source 

term, and is obtained from insertion into A^ LU of the multipole expansions of h and its space- 
time derivatives. Note that M.(T pv ) = since T pv has a compact support. We want to 
compute the multipole moments of an extended post- Newtonian source (one for which the 
PN approximation is physically meaningful). To this end we first consider the following 
quantity: 

r B M(A^)] , (2.5) 



A"" = h" u - FP D^ 1 

B=0 



which is made of the difference between W, the solution of the field equations ()2.2|) valid 
everywhere inside and outside the source, and a particular object obtained from the retarded 
(TV) integral of the multipole (M) expansion of the gravitational source term A pv . Here the 
retarded integral means the usual flat space-time expression 

,-, i ,/ x 1 f d 3 x.' J . |x — x'|\ . 

Since Eq. (J2.6j) extends up over the whole space, x' £ IR 3 , including the region inside the 
source where the multipole-moment expansion is not valid (it diverges at the spatial origin 
|x'| — > located inside the source), one is not allowed to use directly the retarded integral 
as it stands. This is the reason for the introduction in the second term of Eq. ()2.5|) of a 
particular regularization process defined by the finite part when a complex number B tends 
to zero (this operation is abbreviated as FP# = o), and involving the regularization factor 



r 



B 



^|x|»^-j , (2.7) 

which is to be inserted in front of the multipolar-expanded source term of the retarded 
integral. Here ro denotes an arbitrary constant scale having the dimension of a length. 
Since the divergence of the retarded integral is at the origin of the coordinates, |x'| — > in 
(12. 6|) . the constant r plays the role of an ultra-violet (UV) cut-off in the second term of Eq. 
(12. 5|) . However we shall see that in the expression of the multipole moments themselves, 
given by Eqs. (j2.11j) or (j2.13J) below, the same constant ro will appear to represent in fact 
an infra-red (IR) cut-off. 

From Eq. (|2.5|) . and noticing that the second term is already in the form of a multipole 
expansion, we can write the complete multipole decomposition of the external field as 

M{h" u ) = M{W) + FP D^ 1 k B M{A^)] . (2.8) 

B=0 L 
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This is a solution of the vacuum Einstein field equations, 



UM(hr) = M^), (2.9a) 
d v M{hT) = 0, (2.9b) 

valid in the exterior of the source where M.{T^ U ) = 0. The first term in Eq. (|2.8J) is given 
by an homogeneous solution of the wave equation: CLM(A MI/ ) = 0. The second term in (J2.8|) 
represents a particular, inhomogeneous, solution, which arises because of the nonlinearities 
in the external gravitational field, and can be computed by means of the multipolar-post- 
Minkowskian algorithm of Ref . j40j . In the present paper we shall not consider the second 
term in ()2.8|) because its contribution, encompassing many non-linear effects, has already 
been computed in jijj up to 3.5PN order for compact binaries. We shall define our source 
multipole moments from the contribution _M(A M ^), which can be viewed in fact as the 
"linearized" part of the multipolar decomposition ([2.8)1 . 

It has been proved in Refs. 37, |3^, 3^] that: (i) the multipole expansion A4(A 11U ) of 



the quantity defined bv Eq. (12. 5 | ) can be computed using the standard formulas (given 
for instance in Refs. 4j| and |43| ) for the multipole expansion outside a compact support 



source; (ii) the multipole moments admit a very simple expression in the case where the 
matter source is slowly-moving (existence of a small PN parameter e ~ v/c). The result we 
find reads as 4 



AC + °° (-Y 

1=0 



-W v (t-- 



r ' L \ 



(2.10) 



where the time- dependent functionals TO^ so introduced, which depend on the retarded 
time u — t — r/c, are explicitly given by 5 



; M ) = FPy d 3 x|xJ B x L r^(x, M ). (2.11) 



The integrand of the multipolar functional (j2.11j) involves the post-Newtonian expansion 
of the total pseudo stress-energy tensor given by (|2.3jl . where the formal operation of taking 
the PN expansion is denoted by means of an overbar, i.e. T^ u = PN [t^ u ]. This is the crucial 
point on which we recognize that the expression ()2.11|) is valid only for extended PN sources, 
whose compact support extends well within their own near zone (see |38| for details). The 
other important feature of Eq. (|2.11|) is the presence of the finite part operation when 
B — > 0, with regularization factor given by ()2.7|) . The role of the finite part is to deal with 
the IR divergences initially introduced in the multipole moments by the fact that the PN- 
expanded integrand of the multipole moments diverges at spatial infinity (when r — > +oo). 
By contrast, we recall that the finite part in the second term of (12.5)1 was to take care of 
the UV divergences when r — > 0. The fact that the same finite part operation appears to be 
either IR or UV depending on the formula is made possible by the properties of the complex 



4 The notation is: L = i\ ■ ■ ■ ii for a multi-index composed of £ multipolar indices i\, • • • , it, dt = di 1 ■ ■ ■ di e 
for the product of £ partial derivatives di = d jdx l \ similarly xl = x% x ■ ■ • %i e for the product of £ spatial 
vectors x l = xi. In the case of summed-up (dummy) multi-indices L, we do not write the I summation 
symbols, from 1 to 3, over their indices. 

5 With a slight abuse of notation the generic source-point on which one integrates in p.llfl is denoted by 
x which is not the same as the field-point appearing in the R.H.S. of Eq. I|2.10|l . 



7 



analytic continuation (with respect to B e C). We also mention the fact that the two terms 
in Eq. ()2.8|) depend separately on the length scale r , but that this dependence is in fact 
fictitious because the ?Vs can be shown to cancel out. [To prove this the best way is to 
formally differentiate the R.H.S. of Eq. (|2.8j) with respect to r .] 



B. The STF source multipole moments 

In the present approach it is convenient to work with the equivalent of the multipole 
expansion (j2.10|) - (|2.11|) but written in symmetric and trace-free (STF) guise. We present 
only the results. For the multipole expansion we have 



dr 



-JT(t- r - 



(2.12) 



where the multipole moment functionals Tj^{u) are now STF with respect to their I indices 
L = i\ ■ ■ ■ ie. The JF^'s differ from their counterparts Ti!£ parametrizing the non-STF 
multipole decomposition (j2.10j) . They are given by 6 



FP J d 3 x\3i\ B x L j dz8i(z)T tlv (y,u + z\x\/c). (2.13) 



Equation ()2.13j) involves an extra integration, with respect to its non-STF counterpart (|2.11j) . 
over the variable z, and with associated "weighting" function 8^(z) given by 



5Az) = 



dz 6f(z) 



lim 8f(z) 



(21+1) 
1, 

5{z). 



'1-z 



2\i 



(2.14a) 

(2.14b) 
(2.14c) 



Here 5(z) is Dirac's one- dimensional delta-function. 

To obtain the source multipole moments, we decompose the function T*£ ', which is al- 
ready STF in its I indices composing L, into STF irreducible pieces with respect to all its 
spatial indices, including those coming from the space-time indices \w = {00,0i,zj}. The 
appropriate decompositions read 7 



Rj 



T 0i _ (+) T . , . . (0)t i , x. .(-)Tt 1 



> i 



^U ijL + STF STF W +1) f/ aiL -i + 6 lH ^a 



jL-1 



+ 5nV L 



(2.15a) 
(2.15b) 



(2.15c) 



6 We denote the symmetric-trace-free (STF) projection by means of a hat, xl = STF(xi 1 
times by means of brackets () surrounding the indices, xl = 

7 We denote by £yfc the usual Levi-Civita anti-symmetric symbol such that £123 = +1. 
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We have introduced ten STF tensors Rl, ^Tl+i, • • ■ , (~ 2 >Ul-2, Vl, equivalent to the ten 
components of the original tensor ' . Because of the harmonic-gauge condition (j2.9b|) . only 
six of these tensors are independent, and we are led to a set of six source multipole moments, 
denoted {II, Jl, Wl, X l , Y l , Z l }. These moments are defined in such a way |38| that the 
four last ones, {Wl, Xl, Yl, Zl}, parametrize a mere linearized gauge transformation of the 
"linearized" part of the multipolar metric, and consequently do not play a very important 
role. In practice the moments {Wl, ■ ■ ■ , Zl} appear only at high PN order, where they 
can be typically computed with Newtonian precision, so they do not pose computational 
problems. They have already been taken care of in paper I. 

The "main" source multipole moments are the mass-type moment II and the current- 
type one Jl- In Section ITTTI we shall concentrate our attention on the mass moment II with 
full 3PN accuracy. Having in hands the STF-irreducible decompositions (|2.15jh we obtain 
in the generic case where £ > 2 (i.e. for non-conserved, arbitrary time-varying, moments): 

1 _ * 4 



h = ^ + 3V t ) - Hii + WTWTJ) ™U L , (2.16a) 

(2.16b) 



Jl ic L+ 2£c 2 



where time derivatives are indicated by dots. To express the results (j2.16J) in the best way, 
we introduce the following definitions, 



7=00 _|_ jdi 



C 2 



(where r** = 8^) , (2.17a) 



= — , (2.17b) 



c 



r i] . (2.17c) 



For simplicity's sake we omit the overbar of the E^'s to indicate the post-Newtonian ex- 
pansion, but we do not forget that these quantities are given by, and should be treated as, 
formal PN-expanded expressions. The STF source moments, for multipolarities £ > 2, are 
then given by [iH 

i l( u) = fp / a ixi» /_' «b [ Sl iL S - c2( /f j+a. 3) S l+1 iiL i, 

2(2£ + l) ■■ 1 

+ c*(£ + l)(£ + 2)(2£ + 5) 5w ^ L ^S ^ u + z ^ c ^ ( 218a ) 



B=0 



J L {u) = FP e ab<ie I d 3 x|x| B / dzl 5ex L -i >a 



i 

-l 



2£ + 1 • 

■5t+iXL-i>acEbc K x ' u + z\x\/c) , (2.18b) 



c 2 (£ + 2)(2£ + 3) 



where the E^^'s are evaluated at the position x and at time u + z\x\/c. In the limiting 
case of linearized gravity, we can replace t^ u by the compact-support matter tensor and 
ignore the finite part procedure (FPb=o), so we recover the linearized-gravity expressions 
obtained in Let us emphasize that Eqs. ()2.18|) are u exact" , in the sense that they are 
formally valid up to any PN order. In practice, the PN-expanded moments ()2.18|) are to be 
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computed by means of the infinite post-Newtonian series 
J dz Sg(z) £(x, u + z|x|/c) 

Oik,£ 

In a separate work j^H we shall derive some alternative expressions of the PN moments 
fl2.18j) in the form of integrals depending only on the boundary at infinity {i.e. |x| — > +00, 
u = const). 



x| d 

c du 



+00 / |__| Q N 2k 

k=0 

(21 + 1)!! 
{2k)\\{2e + 2k + 1)!! 



£(x, m) 



(2.19a) 
(2.19b) 



C. The conserved monopole and dipole moments 

In the case of the non-radiative moments, i.e. the mass monopole M (£ — 0) and the 
mass and current dipoles Mj and Si (£ = 1), things are a little bit more involved than what 
is given by Eqs. (|2.18|) . The monopole and dipoles are conserved by virtue of the source's 
equation of motion, Eq. ()2.4|) . namely 

M = 0, (2.20a) 
Mi = 0, (2.20b) 
Si = 0. (2.20c) 

In particular M denotes the ADM mass of the source. As shown in [38| the conserved 
monopole and dipoles can be written into the form 

M = I + 51, (2.21a) 
Mi = Ii + SIi, (2.21b) 
St = J t + 5Ji, (2.21c) 

where the first pieces J, Jj and Ji are defined by the same formulas as Eqs. ()2.18|) but in 
which we set either i = or I = 1, and where the extra pieces follow from Eqs. (5.6) in [38], 
together with (5.4) and (4.5) there, and are explicitly given by 

51 = FP B J d 3 x|xf^ J + ^8 x x h Y,A , (2.22a) 

6Ii = FP B j d 3 x|xf ^_ J ^dz ^XiE^ - S ^ 2) + ^5 2 x ib E ab \ , (2.22b) 

5 J, — FP B J d^\k\ B e iab ^ j dzS^-V. (2.22c) 

Time anti-derivatives are denoted by superscripts (— n); Sq and Si refer to the function given 
by ()2.14j) : like in Eqs. ()2.18)1 the integrands are evaluated at point x and at time u + z|x|/c. 
The quantities 51, 5Ii and 5Ji are precisely such that the "total" moments M, Mj and Si 
obey the conservation laws ()2.20|) as a consequence of the matter equations of motion (see 
Ref. J38[ for further discussion). 
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The chief feature of the expressions (|2.22jl is that they involve an explicit factor B in 
front, and therefore they depend only on the behavior of the integrand when |x| — > +00, 
since they are zero unless the integral develops a pole ~ 1/B due to the behavior of the 
integrand near the boundary at infinity. We shall give more details on the way we compute 
such integrals "at infinity" in Section II V Dl 



III. THE MASS MULTIPOLE MOMENTS AT THE 3PN ORDER 



In this Section we derive the mass-type source multipole moment Ii (for arbitrary £ > 2) 
at the 3PN approximation, for general extended PN sources. From Eqs. (J2.18|) we see that 
one must obtain E with full 3PN accuracy, but E, at the 2PN order only, and at 1PN 
order. For this purpose, we make explicit the components of the E^'s, defined by Eqs. 
()2.17|) . in terms of a certain set of retarded "elementary" potentials: V, Vi, Wij, X, Ri, Zy, 
solutions of appropriate iterated d'Alembertian equations. Although devoid of any direct 
physical meaning, these potentials have proved to constitute some very useful "building 
blocks" for practical PN calculations on grav itational- wave generation (paper I), as well 



. In paper I we systematically expanded 

potentials 



as in the problem of equations of motion 2 

all the retardations in V, Vj, • • • and introduced some associated "Poisson-like 
U, Ui, • • • . Here we shall come back to the same retarded-like potentials V, Vi, W^, ■ ■ ■ , Zy 
as in the equations of motion; of course we are motivated by the fact that they have already 
been computed in So we shall redo entirely the computation of paper I, using different 
elementary potentials and also more systematic Mathematica programs, and in the case of 
general orbits. Our results will match perfectly with those of paper I. 

Let us denote the "matter" parts in the total density contributions (j2.17|) by 



a. 
a 



^00 _|_ J^l 



a 



I = SijT 3 



T 



0/ 



c 

rpij 



(3.1a) 

(3.1b) 
(3.1c) 



Our chosen definitions for the elementary retarded-type potentials, which involve non-linear 
couplings appropriate to 3PN order, are 



V 



w, 



X 



Ri 



Zij 



D^ 1 [-47rGa] , 

AnG(a ij - 5 i: j(j kk ) - diVdjV 



□- 1 



□- 1 



-AtcGcjuV + Wijd^V + 2VidAV 
+ Vd 2 t V + ^(d t V) 2 - 2d i V j d i V i 



-AnG^V - aVi) - 2d k VdiV k - -d t VdiV 

-AnG(a tJ - 5 i0 a kk )V - 2d {j yd t V 3) 
+ diVtdjVk + d k Vd k V 3 - 28^8^ 



(3.2a) 
(3.2b) 

(3.2c) 

(3.2d) 
(3.2e) 
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- 5 tJ d k V m (d k V m - d m V k ) - -^{dtVf^ , (3.2f) 

together with the spatial traces denoted by W = Wu and Z — Notice that we shall not 
need some higher-order potentials, called T and Y, which were crucial in the 3PN equations 



of motion |22|. The 3PN multipole moments are in this sense "less non-linear" than the 
3PN equations of motion. 

Like in paper I we find it convenient to decompose II into three pieces corresponding 
to the three terms in ()2.18af . respectively referred to as "scalar" (S), "vector" (V) and 
"tensor" (T). Applying the formula (j2.19j) we further decompose each of these pieces into 
parts, labelled as I, II, III and so on, according to the successive PN contributions. Hence, 
we write 

I L = SI L + SII L + SIII L + SIV L 
+ VI L + VII L + VIII L 

+ TI L + TII L , (3.3) 

in which we consistently neglect all terms that are higher-order than 3PN. 8 Without further 
comment and proof, we give the explicit expressions of all these separate pieces, which are 
equivalent to the similar expressions given by Eq. (4.2) in paper I. Concerning the "S-type", 



SI; 



FP / d 3 x|x| B x L < a - 

B=0 / 



ttGc 4 

4 « 1 
+-t^d i V i d t W ii - -^-Z^V - -4-dAVRi 



vrGc 6 
8 



+^rrd i V j d j R i - 



2 



sn. 



ttGc 6 
1 

1 



ttGc 6 



— 6 A(VZ) ,0.4a) 



d 2 



— — r FP — / d 3 X \k\ B { |x| 2 X L 

2c 2 (2£ + 3) b=o dt 2 J 



AV 



ttGc 4 



VidtdiV 



ttGc 4 



WijdijV 



1 



2nGc 4 



{d t Vf + 



We generally do not indicate the PN remainder term O (c 7 ) 
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2 -^*lV 2 - ^ 2 d t [d^V 2 )^ - V 2 d t (\x\ 2 x L ) 
2£ + 3 



^^d t [di{VW)\x\ 2 x L - VWd t (\x\ 2 . 



■x L ) 



SIII Z 



SIV L = 



21 + 3 - 1 

\ FP 

8c 4 (2£ + 3)(2£ + 5) b=o 

2(2£ + 5). .„. __n 



dt* 



cfx \x\ < |x| 4 x L a 



:(|x| 4 £l)' 



^ + 5 ^ lxl 2 £ f \/ 2 - d \d (V 2 )\x\ 4 x r - V 2 d- n - 1 

1 F p f_ [ d 3. 

48c 6 (2£ + 3) (2£ + 5)(2£ + 7) s=o d*e 7 



! x |x| s |x| 6 fL(7 . 



Then, the vectorial V-parts are 



VI L = 



4 ^ +1 ) pp 1 f A 
c 2 (£+l)(2£ + 3) b=o dt 7 

1 3 



:\B 



(2 2 



4 



a.V 2 



12 2-2 1 

c 4 c 4 JJ c 4 J 3 c A 3 3 2nGc A 1 



VII L = 



^ A ^) " 2^? A ^> + 2^ A( ^} ' 



4 (£ + 



2(2^+1) Fp 
■ l)(2£ + 3)(2£ + 5) b=o 
2 2 1 



^1 
dt 3 

+ ■ 



d 3 x Ixl 



|x| 2 x iL 



2£ 4- 5 1 r 



yjTT . . 2f+l 

L ~ 2c 6 (£+l)(2£ + 3)(2£ 

Finally the tensor parts read 



! X,L- V^(|x| 2 X iL )] j 

FP ^ / d 3 x|xf:r lL |x| 4 o 

+ 5)(2£ + 7) s=o dt 5 J 11 ^' 1 



TI L = 



2 ^ +1 ) fp^. /> X |X, 

^ + l)(£ + 2)(2£ + 5) b=o d* 2 y 1 



X 



ij 



2 criVj 
c 2 
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,k md > v * + d? d - VAV ' - sb™} • (3 ' 6a) 



THr = 2i + 1 FP — f d 3 xlxJ B x-rlxl 2 <!V- 

c6(£+l)(£ + 2)(2£ + 5)(2£ + 7)Bio d t 4 y dx l x l x ^l x l \^ 

+ 4^G^ y }- (3 - 6b) 

These formulae stricto sensu are valid for general time-varying multipole moments having 
£ > 2. However they constitute also the main contributions in the conserved monopole and 
dipole moments [I = 0, 1) as well. In fact we shall prove in Section IV Bl that Eqs. (|3.4jl — (|3.6|) 
already give the correct answer for the 3PN mass dipole moment of point particle binaries, 
i.e. Mi = Ii and the quantity Sli given in Eq. (j2.22b|) is zero at 3PN order. 



IV. HADAMARD REGULARIZATION OF THE MULTIPOLE MOMENTS 

We now specialize the general expression of the 3PN mass moments to compact binary 
systems modelled by point particles. To this end the first task is to compute all the necessary 
potentials {V, V{, Wy, ■ • •}, in the case of delta-function singularities, using Hadamard's 
regularization. Actually the computation of all these potentials has already been done at 



the occasion of the 3PN equations of motion, and we refer to [22( for the details. Our next 
task is to insert these potentials, and their space-time derivatives, into Eqs. ()3.3p - (j3.6j) for 
the quadrupole and dipole moments, following the prescriptions of the Hadamard or more 
precisely the pure-Hadamard-Schwartz (pHS) regularization. 

A. Pure Hadamard Schwartz regularization 

Let us first recall the two concepts that constitute the basis of the "ordinary" Hadamard 
regularization (13 . 141 . 9 The first one concerns the partie finie of a singular function at 



the value of a singular point. The generic function we have to deal with reads -F(x), where 
x G M 3 , and becomes singular at the two point-particle singularities located at the positions 
yi and y2 (in the harmonic coordinate system). The function -F(x) is smooth (C°°) except 
at yi and y2, and admits around these singularities some Laurent-type expansions in powers 
of ri = |x — yi\ or r 2 = |x — y a | . When r\ — > we have, ViV G N, 

H*)= E ^/ P (ni) + o(rf), (4.1) 

where the Landau o-symbol takes its usual meaning, and the coefficients 1 / p (rii) are functions 
of the unit vector ri! = (x — yi)/ri. 10 We have p G Z, bounded from below by some typically 
negative inte ger p n depending on the F in question. The class of functions such as F is called 
T\ see Ref. |23( for a fuller account of the properties of functions in this class. Now the 



^^fer to for a digest of possible variants of Hadamard's regularizatioon. 

10 For clearer reading, we use a left-side label 1 like in i/ p (ni) when the quantity appears within the text, 
however the label is always put underneath the quantity when it appears in an equation like in H4.ll) . 
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Hadamard partie finie of F at the singular point yi, denoted (F)i, is defined by the angular 
average 

(4-2) 

where dQ\ = dfl(ni) is the solid angle element on the unit sphere centered on y^. Note 
that the spherical average (|4.2j) is performed in a global inertial frame. In the context 
of the extended-Hadamard regularization j24| one defines the regularization ()4.2jl in the 
Minkowskian "rest frame" of each particle. A distinctive feature of the partie finie (|4.2|) is 
its "non-distributivity" in the sense that 

(FG)i ^ (F) 1 (G) 1 in general for F,G G T. (4.3) 

The second notion in Hadamard's regularization is that of the partie finie of a divergent 
integral, which attributes a value to the integral over M 3 of the function F(x). Consider first 
two "regularization volumes" around the two singularities yi and y2. We can specifically 
choose two spherical balls (in the considered coordinate system), Bi(s) and B 2 (s), centered 
on the singularities, each of them with radius s. The integral of F over the domain exterior 
to these balls, i.e. R 3 \ Bi(s) U B 2 (s), is well-defined for any s > 0. Hadamard's partie finie 
(Pf) of the generally divergent integral of F is then defined by the always existing limit 

Pf sliS2 / d 3 xF(x) = limi / d 3 xF(x) 



R 3 \Bi(s)UB 2 (s) 



^^(S^KS^ 1 - 2 }- (4 - 4) 

The extra terms, which involve some parties finies in the sense of (J4.2j) . are such that 
they cancel out the singular part of the "exterior" integral when s — ■> 0. Here the symbol 
1 «-> 2 means the same terms but corresponding to the other particle. The two constants s\ 
and s 2 entering the logarithmic terms of this definition play a very important role at 3PN 
order. A way to interpret them is to say that they reflect the arbitrariness in the choice 
of the regularization volumes surrounding the particles. Indeed it can be checked that the 
Hadamard partie finie (|4.4jl does not depend, modulo changing the values of si and s?, on 
the shape of B\ and B 2 , above chosen as simple spherical balls (see the discussion in |23|]). 

The two notions of partie finie, (|4.2jl and (|4.4jl . are intimately related. Notably the partie- 
finie integral (|4.4j) of a gradient is in general non-zero but given by the partie finie, in the 
sense of (J4.2|) . of some singular function (see 23|| for more details). With the definitions (|4.2|) 



and (|4.4j) one can show that, if we want to dispose of a local meaning (at any field point 
x) for the product of F with a delta-function, say F(x) 5(x — yi), then one cannot simply 
replace F in front of the delta-function by its regularized value. This is a consequence of 
the non-distributivity of Hadamard's partie finie, Eq. (|4.3J1 . Thus, 

F(x)5(x-yi) ^ (F)i<5(x-yi) in general for F E T. (4.5) 

It is quite evident that the two properties (J4.3)) and (J4.5|) are problematic. A remarkable 
fact is that the problem of the non-distributivity, Eqs. (|4.3|) and ([4.5)1 . arises precisely at 
the 3PN order, for both the radiation field and the equations of motion, and not before that 
order. In the problem of the equations of motion we could deal with the properties ()4.3|) and 
(14.5)1 by implementing the extended Hadamard regularization of Refs. [23, 24]|. We have 
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not (yet) succeeded in applying the extended Hadamard regularization to the problem of 
gravitational wave generation. For the present paper we choose to follow a different route, 
and adopt the pure- Hadamard- Schwartz (pHS) regularization defined in Ref. |27j ]. 

The pHS regularization is a specific "minimal" variant of the Hadamard regularization, 
which is designed in such a way that it avoids, by its very definition, the problematic con- 
sequences (|4.3|) and (J4.5|) of the "ordinary" Hadamard regularization. It applies to the case 
relevant here where the singular function, say Fl G T, is made of sums of products of 
the non-linear potentials V , Vi, W^, ■ ■ ■ and their space-time derivatives diV, ■ ■ ■ , and is 
multiplied by some (regular) multipolar factor xl, that is 

F L (x) = x L V[V, V, Wij,---, diV,---}, (4.6) 

where V denotes a certain multilinear form, i.e. a polynomial in each of its variables 
V, Vi, ■ ■ ■ . The rules of the pHS regularization are: (i) an integral J d 3 xFi(x), where 
Fl takes the form (J4.6|) . is regularized according to the partie finie prescription (|4.4|) like in 
the ordinary Hadamard regularization; (ii) we add the contribution of distributional terms 
coming from the derivatives of potentials, d^V, ■ ■ ■ , according to the usual Schwartz distribu- 
tion theory 14] (this point is detailed in Section riVB|) : (iii) the regularization of a product 
of potentials V, Vi, W^, ■ ■ ■ (and their gradients) at a singular point is assumed to be "dis- 
tributive", which means that the value of Fl at point 1 (say) is given by the replacement 
rule 

(f l )i — > tfv[(v)i, ryOi, (w^, ■■■ , -••], (4.7) 

where the partie finie ()4.2|) is applied individually on each of the potentials, and = 
STF^ 1 • • • yY); and (iv) a "contact" term, i.e. of the form F L (x.) <5(x — yi), appearing 
in the calculation of the sources of the non-linear potentials and corresponding to their 
"compact-support" parts, is regularized by means of the rule 

F L (x)<5(x- yi ) — y^V[{V) x , (VJ)i, (Wi;)i, • • •] 5(x - yi) . (4.8) 

The rules (14. 7j) and (|4.8|) of the pHS regularization are well defined, and are not submitted, by 
definition, to the unwanted consequences of the non-distributivity of the ordinary Hadamard 
regularization: (J4.3)) and (J4.5|) . However, as we shall emphasize in Section |3 the pHS 
regularization becomes physically incomplete at the 3PN order, in the sense that it must be 
augmented by certain ambiguous contributions, which a priori cannot be determined within 
this regularization scheme. 

Our motivation for introducing the pHS regularization is that it constitutes in some sense 
the core of both the Hadamard and dimensional regularizations (27l.l32^. By "core" we mean 
that it will yield the complete and correct result for all the terms but for a few, and for those 
which cannot be determined unambiguously the undetermined part will take in general a 
very special and limited type of structure. Hence the correct result is obtained by adding 
to the pHS result a limited number of "ambiguous" terms, parametrized by some arbitrary 
numerical coefficients called ambiguity parameters. 

In dimensional regularization the undetermined terms correspond exactly to the contri- 
bution of poles oc 1/e, where d = 3 + e is the dimension of space. The complete result in 
dimensional regularization appears therefore as the sum of the pHS result and what we call 
the "difference", namely the pole part ~ 1/e which can be quite easily obtained from the 



De quit' 

expansion near the singularities of the functions involved (27 . l32t] . and which is nothing but 
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the difference between the dimensional and the pHS regularizations. The method for deter- 
mining the ambiguity parameters is therefore to equate the ambiguous terms, as they are 
defined with respect to the pHS regularization, to the latter difference. In the present paper, 
we compute the pHS regularization of the 3PN binary's quadrupole moment; this constitutes 
the first and necessary step toward the complete calculation by dimensional regularization 
(see Ref. j32| for a summary of the method). 



B. Schwartz distributional derivatives 



We detail here an important feature of the pHS regularization, namely the systematic 
use of distributional derivatives a la Schwartz 14]|. Recall first that previous work on the 



equations of motion |22j showed that the Schwartz distributional derivatives yield ill- defined 
(formally infinite) terms at the 3PN order in ordinary three-dimensional space. This was a 
motivation for introducing some appropriate generalized versions of distributional derivatives 
in the context of the extended Hadamard regularization [23j]. However, one can show [2^ 
that, by working in a space with d dimensions instead of 3 dimensions, and invoking complex 
analytic continuation in d, the latter ill-defined terms are in fact rigorously zero. The usual 
Schwartz distributional derivatives are therefore well-defined in the context of dimensional 
regularization, and they have been shown to contribute in an essential way to the final 
equations of motion at 3PN order |27| . 

In the present paper we include the Schwartz distributional derivatives as part of the 
calculation based on the pHS regularization. However, as we just pointed out the Schwartz 
derivatives yield ill-defined terms in 3 dimensions, so we shall compute them in d spatial 
dimensions, and then take the limit 

e = d-3 -> 0. (4.9) 

This permits us to cancel out (by dimensional continuation) all the formally divergent terms 
and to perform a perfectly rigorous calculation. Of course, this way of handling the Schwarz- 
tian distributional terms shows that in fact the calculation will already constitute a part of 
a complete calculation using dimensional regularization. However the spirit is different. In 
the present calculation we use dimensional continuation as a mathematical trick enabling 
us to give a well-defined meaning to a limited number of terms which would be otherwise 
infinite. In a real computation based on dimensional regularization the scope is broader, and 
we should start from the Einstein field equations in d dimensions and consistently perform 
all the derivations for arbitrary d G C before eventually taking the limit (|4.9J1 . In the present 
paper we shall perform our calculation of the Schwartz derivatives based on the expression 
for the multipole moment given by Eqs. ()3.4D - (J3.6I) . i.e. without taking into account the 
modification of the various coefficients which would come from the Einstein field equations in 
d dimensions. The result in the limit e — > will however exactly be the same as by including 
such d-dependent coefficients because the distributional parts of Schwartz derivatives do not 
generate any poles proportional to \ je. 

One may ask why is it possible to choose, in the context of Hadamard's regularization, 
different prescriptions for the distributional derivatives, and nevertheless obtain the same 
physical result at the end ? For instance what would happen if instead of using the Schwartz 
distributional derivative in the way we have just described, we adopt the generalized deriva- 
tives of the extended Hadamard regularization |23[ ? The answer which emerges from our 
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detailed computations is that the difference between the final results obtained by different 
prescriptions takes the form of the "ambiguous" terms, which are given in the case of the 
quadrupole moment by the R.H.S. of Eq. (J5.6|) below. Thus, different calculations are in 
fact equivalent modulo some simple redefinition (or shift) of the values of the ambiguity 
parameters £, k and (. 

The d- dimensional calculation of the Schwartz distributional derivatives in the 3PN mo- 
ments essentially necessitates the same ingredients as in the problem of equations of motion 
(2?}. We introduce some elementary Poisson kernels u% and v%, solving the equations 

A Ul = -47r5 (d) (x- yi ), (4.10a) 
Av 1 = ui, (4.10b) 

where A is Laplace's operator in d dimensions and 5^ is the Dirac delta-function in d di- 
mensions. These kernels play a crucial role in the construction of the ci-dimensional versions 
of the non-linear potentials [27]. They parametrize the compact-support potential V at 1PN 
order; evidently Ui enters the Newtonian part of V while the twice-iterated Poisson kernel 
v 1 is used for the 1PN retardation. The kernels are given by 

u x = ~kr\- d , (4.11a) 

k r 4_d 

^ - fco- (4 ' nb) 

where ri = |x — y x | and k is related to the Eulerian T-function by 

t _ r (¥) 



d-2 
TT 2 



(4.12a) 



limfc = 1. (4.12b) 

(i — >3 

The second partial derivative of ux, and the fourth partial derivative of v±, will contain, be- 
sides an ordinary singular function (or pseudo-function) obtained by performing the deriva- 
tive in an "ordinary" sense, a distributional component proportional to 5^ d \ and given by 

4.77- 

dlM = <9JM| M M(x- yi ), (4.13a) 

I ordinary 

d? jH {vi) = d* kl ( Vl )\ -^^faSu + SnSfl + SuS^SMfr-y!). (4.13b) 

lordinary "I" + ^) V ' 

These expressions can be derived as particular cases of the Gel'fand-Shilov formula [45j|. In 
addition, we can treat the distributional tzme-derivatives in a very simple way from the rule 
d t = —v\ di applicable to the purely distributional part of the derivative. 

The expressions ()4.13|) permit for instance the computation of the distributional derivative 
dfjV at the 1PN level to be inserted into the 1PN source term ~ WijdfjV in the expression of 
II, cf. Eq. (|3.4aj) . Let us emphasize that the previous method of introducing the Schwartz 
distributional derivatives in the pHS formalism, i.e. by means of dimensional continuation 
in d, is probably the only rigorous way to do it. An alternative approach would consist 
of staying in 3 dimensions, and employing the generalized derivative operators defined in 
[23l | (they act on singular functions of the class T instead of smooth "test" functions with 
compact support as in Schwartz's distributional theory). But then the result will differ from 
Schwartz's derivatives by some terms having the structure of the ambiguous terms in the 
R.H.S. of Eq. flEHJ). 
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C. Three-dimensional partie finie integrals 



The main basis of the computation of the 3PN multipole moments of point particles is 
to perform explicitly many three-dimensional non-compact support integrals in the sense of 
the Hadamard partie finie (J4.4j) . In addition to the partie finie we have to take care of the 
finite part process based on analytic continuation in B £ C to deal with the boundary of 
the integrals at infinity Therefore we must compute explicitly many integrals of the type 



which depend a priori on the two UV-type length scales S\ and S2 associated with the 
Hadamard partie finie (|4.4|) . and on the IR-type length scale r introduced into the general 
formalism of Section ITD through the regularization factor |x| B = \x/ro\ B . 

The function -F(x) in Eq. (J4.14j) stands for a non-compact support function which, as 
far as its UV properties are concerned, belongs to the class of singular functions J 7 , i.e. 
admits some expansions of the type (j4.1|) . The IR behavior of F(x), when |x| — > +00, will 
be specified below. -F(x) contains also some multipolar factor such as xl but for simplicity's 
sake we do not indicate here the multi-index L. In the general case -F(x) admits an expression 
such as Eq. (14. 6|) . i.e. it is given by some multi-linear functional of the elementary potentials 
V, Vi, Wij, ■ ■ ■ and their derivatives. The function F represents the sum of all the non- 
compact support terms in Eqs. (|3.4j) - (|3.6|) . taking into account only the ordinary parts of the 
derivatives. The compact support terms in (|3.4|) - (|3.6|) . as well as the purely distributional 
parts of the Schwartz derivatives [calculated with (|4.13|) ]. are treated separately using the 
rules of the pHS regularization for contact terms, see Eqs. (|4.7|) - (j4.8|) . On the other hand, 
we shall point out in Section QV Dl that for many non-compact terms in Eqs. (|3.4|) - (j3.6|) we 
can after integrating by parts perform a much simpler computation of these terms, confined 
to the boundary of the integral "at infinity" and depending on the sole properties of the 
finite part operation FPb =0 . 

In this Section we explain our practical method for dealing with the integral (j4.14|) . The 
basic idea is to relate (j4.14|) to an integral which is convergent at infinity, on which we can 
thus remove the finite part at B — 0, and then to compute this integral by means of the 
very efficient method of "angular integration" described by Eq. (4.17) in Ref. j23{. We 
assume for this calculation (this will always be verified in practice) that F admits a power- 
like expansion at infinity, when r\ — > +00 with t = const, of the following type (for any 
large enough M) 



where the coefficients i<£>& depend on the unit vector rii = (x — yi)/r\. The index k is 
bounded from below by some ko £ Z. For convenience we have singled out the singularity 1, 
and considered the expansion when r\ — > +00, instead of the more natural choice |x| — > +00. 
Introducing such an asymmetry between the points 1 and 2 is only a matter of convenience, 
but in fact it is quite appropriate in the present formalism because we shall later use the 
method of "angular integration" |23| which already particularizes the point 1, around which 
the angular integration is performed. An advantage is that a good check of the calculation 
can be done at the end since the final result will have to be symmetric in the particle 
exchange 1 2. Next we define an auxiliary function iF^ by subtracting from F all the 




(4.14) 




(4.15) 
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terms in its expansion (j4.15|) which yield some divergencies at infinity, i. e. 

Foo(x) = F(x)- "4^( n i)- ( 4 - 16 ) 

1 z — ' rv i 

fc <fc<3 1 



The integral of is easily seen to be convergent at infinity, and therefore it can be 
computed with the ordinary Hadamard partie finie prescription given by (|4.4j) . Inserting 
(ETTKJ) into (gJH) we obtain 

X[ Sl ,s 2 ,r ] =Pf Sl , S2 / rf 3 xFoo+ If { Pfsi ' S2 / d " X l 5 | B ^? fe ( ni ))- ^ 4 - 17 ) 

1 fc <fc<3 L •/ r l J 

We transform the second term in the R.H.S. of ()4.17jl . The integrand in the curly brackets 
is replaced by the following equivalent when B — > limited at first order in B, 



J(Px\x\ B ^<p k (n 1 ) = Jd 3 ri r B 1 + B In (^j + O (B 2 ) 



\<Pk(*i). (4-18) 



(Notice our change of integration variable, from x to ri = x — yi; we pose r\ = Ti/tq.) 
This will turn out to be sufficient for our purpose because at 3PN order one can show that 
there are no multiple poles in B, therefore the neglected terms O (B 2 ) will never contribute 
at this order. We can prove that the first term in the R.H.S. of (j4.18j) is zero except when 
k = 3, in which case the integral admits a pole, and its finite part depends on the logarithm 
of the ratio between ro and s\. Furthermore, the next term, carrying an explicit factor B, 
is found to be zero in the case k = 3, so we get 

Fp|pf SliS2 y d 3 x |xf I^ 3 ( ni )j = ln (^) / ^i<f 3 (ni). (4.19) 

Consider next the generic cases where k < 2. It is clear that the integrals are then convergent 
when ri — > so we can ignore the Hadamard partie finie (Pf). By analytic continuation 
in B we find that the first term in (J4.18)) is now zero, and it remains the next one, which 
can clearly contribute only if the integral develops a simple pole ~ 1/B at infinity. When 
t*i — > +oo we have the expansion 



m^He 2 ^. (420) 

m=l 1 



''1 



where the various coefficients a m depend on and also on y 1; and are related to the 
Gegenbauer polynomial C^(t) by 11 



a m (ni) = - yi 



d 

dfi 



(giyQ 
lyil 



(4.21) 

fi=0 



where one sets \x = after differentiation of C^(t) with respect to its argument /i. One may 
want to express (|4.21|) in a more detailed way with the help of Rodrigues' formula for the 



ii 



Here we follow the standard convention for the Gegenbauer polynomial In Eq. I|4.21|) it is calculated 
|yi 



at the value t = — (" iyi ^ ; where (niyi) denotes the usual scalar product and |yi| the usual norm 
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Gegenbauer polynomial. Substituting the expansion (|4.2UJ1 - (|4.21|) into Eq. (J4.18)) we are 
finally in a position to obtain the looked-for result 

I[si,s 2 ,r ] = Pf SljS2 J rf 3 x Foo 

+ ln (^) J <Mi fsKi-i^y dnia ro (m) ^ 3 _ m (m). (4.22) 

This formula is systematically employed in our algebraic computer programs. [Of course, 
the sum in the last term is in fact finite because 3 — m > ko] see Eq. (|4.15|) .] As we 
said the first term (partie finie integral) is computed by means of an "angular integration" 
around the particle 1 following the procedure defined by (4.17) in Ref. [23j. We have found 
that the 3PN quadrupole moment resulting from the systematic application of Eq. (|4.22|) 
is in perfect agreement with the result of paper I, which was derived by "case-by-case" 
integration, i.e. using different methods depending on the type and structure of the various 
terms encountered in the problem. 



D. Contributions depending on the boundary at infinity 

The result (I4.22|) can be applied to any of the non-compact support terms in ()3.4)) - (J3.6)) . 
However, we now show that many terms can be re-expressed, after suitable integration by 
parts, in the form of a surface integral at infinity r = |x| — > +oo. Evaluating the surface 
integral is in general much simpler than performing the "bulk" calculation following Eq. 
(14.22)1 . The first type of term in (|3.4)) - ([3.6j) for which a computation "at infinity" is possible 
takes the form of the finite part (FPb=o) of an integral involving the product of a multipolar 
STF factor x L with the Laplacian of some G G J 7 , having the structure of a product of 
elementary potentials, i.e. 

J L = FP J rf 3 x |x| B x L AG . (4.23) 

There are many such terms in ()3. 4)1 - 1)3.6)) . see for instance the last six terms in Eq. (j3.4aj) . 
The second type of term which is amenable to a treatment at infinity is composed of the 
divergence of some vectorial function Hi G J 7 , containing itself some multipolar factor xl 
(not indicated in our notation for Hi), say 

JC = FP J d 3 x |x| B diHi . (4.24) 

An example is given by the last two terms in Eq. (|3.4b|) . We deal with these two categories 
of terms, Jl and /C, in turn. 

The Laplacian in Jl is integrated by parts, and we are allowed to cancel out the all- 
integrated term which is zero by analytic continuation in B G C (because it is zero in the 
case where dt(B) is chosen to be a large enough negative number), thereby obtaining 

J L = FP [ d 3 x A (|x| B xl) G = FP B(B + 2£ + 1) [ d 3 x |x| B ~ 2 — ^ G . (4.25) 

B=0 J B=0 J r 
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The presence of the factor B means that the result depends only on the polar part ~ 1/B 
of the integral at the boundary at infinity. Since the pole comes exclusively from a radial 
integral of the type J drr B ~ x = r B /B, we need only to look for the term of the order of 
r -t-i j n expansion of G when r — > +00. We compute the expansion of G and obtain 

G = --- + ^-X e (n) + o(-^-] , (4.26) 



where the dots indicate the terms having different magnitudes in 1/r and which thus do not 
concern us for the present calculation. The interesting coefficient in ()4.26|) is JJQ(n), and in 
terms of it we get 



J L = FP B{B + 2£ + 1) r B / dr r B ~ x / dtt h L X e {n) , (4.27) 
B=0 Jn J 

in which we indicated that the radial integral depends only on a neighborhood of infinity, 
from some arbitrary radius TZ up to +00. This point is actually not completely obvious at 
this stage and must be justified in the following way. In the general formalism of Ref. [38], 
which is valid for extended smooth matter distributions, any integral having a factor B in 
front will depend only the behavior of the integrand at infinity. Indeed since the matter 
source is smooth the near zone part of the integral is convergent, thus no poles ocl/B can 
arise due to the UV behavior of the integrand and only the IR-type poles can contribute 
to the value of the integral. When applying the formalism to point particles one must keep 
this feature in mind, and replace the stress-energy tensor of an extended source by the T^ u 
of point-particles in the term in question already in the form, by the previous argument, of 
some far- zone integral. Thus, even for point particles the term depends only on the boundary 
at infinity and does not explicitly involve UV-type divergencies, although it may implicitly 
contain some contributions coming from UV divergencies occuring at previous PN iteration 
steps. Finally, from Eq. ()4.27p we readily find the result 

J L = -{2£ + l) J dflh L X e (n). (4.28) 

We notice that this result is independent of the arbitrary scale TZ introduced in (|4.27j) . as 
well as of the IR constant r . Concerning the integral K defined by Eq. (|4.24|) we proceed 
similarly by integration by parts. We find that the term depends only on the part in the 
expansion of Hi at infinity which goes like 1/r 2 , hence we look for the coefficient ^(n) in 

H l = --- + ±- 2 Y l (n) + o(±) , (4.29) 



and we obtain the simple result (independent of TZ and r ) 

JC = j dttriiYiin). (4.30) 

In summary, many non-compact support terms in (|3.4|) - (j3.6|) . having the structure of Jl 
and /C, are computed by surface integrals at infinity using the properties of the analytic 
continuation in B. The only task is to look for the relevant coefficients in the expansions 
of the integrands at infinity, ()4.26|) or ()4.29|) . and to perform the surface integrals ()4.28|) or 
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([4.30)1 . This saves a lot of calculations with respect to the "bulk" calculation of the Hadamard 
partie finie based on the form found in Eq. (|4.22|) . Of course, the two calculations, at infinity 
and in the bulk, will completely agree, but notice that for this agreement to work, one must 
crucially take into account in the bulk calculation, in addition to the formula (|4.22|) . the 
contribution of the distributional part of derivatives. Thus the Laplacian in Eq. (|4.23|) is to 
be considered in a distributional sense. 

V. MULTIPOLE MOMENTS OF POINT PARTICLE BINARIES 
A. The 3PN mass quadrupole moment 

We have computed the 3PN mass quadrupole moment of point particle binaries, for 
general orbits, using the expressions ()3.4|l - (|3.6jl with t = 2, following the rules of the pHS 
regularization, notably the way ()4.8|) one handles the compact-support "contact" terms, 
and the techniques reviewed in Sections IIV CI and IIV Dl to compute three-dimensional non- 
compact-support integrals. We denote by J? HS [si, s 2 , r ] the result of such pHS calculation, 
in which Si and s 2 denote the two UV cut-offs and r the IR length scale. These constants 
result from the computation of non-compact-support integrals and are shown in our basic 
formula (|4.22|) for the Hadamard partie-finie integral. Now it was argued in paper I that 
the Hadamard regularization of the 3PN quadrupole moment is incomplete, and must be 
augmented, in order not to be incorrect, by some unknown, ambiguous, contributions. 

The first source of ambiguity is the "kinetic" one, linked to the inability of the Hadamard 
regularization to ensure the global Poincare invariance of the formalism (we are speaking here 
of the ordinary or pHS variants of the Hadamard regularization, as well as of the "hybrid" 
regularization which has been used in paper I for the generation problem 12 ). As discussed 
in Section X of paper I we must account for the kinetic ambiguity by adding "by hands" 
a specific ambiguity term, depending on a single ambiguity parameter called (. Following 
here exactly the same reasoning we add to the pHS result the same type of ambiguous term, 
which means that we must consider as correct the following 3PN quadrupole moment, 

I ij [s 1 ,s 2 ,r ;(] = /^ HS [si,s 2 ,r ] 

+ -±— ^^f + 1^2, (5.1) 

where the extra term, purely of 3PN order, involves an unknown coefficient £. Here the coor- 
dinate velocity is denoted v\, and the factor 44/3 is chosen for convenience. The parameter 
( will turn out to be different from the parameter ( of paper I because we are adding it to 
the result of the pHS regularization, instead of the "hybrid" Hadamard-type regularization 
considered in paper I. [The hybrid regularization differs from the pHS one by the way the 
contact terms are computed, which takes into account the properties of non-distributivity 
(|4.3|) and (|4.5|) and is more like the one of the extended Hadamard regularization [23j], and 



An exception is the extended Hadamard regularization which is in principle able to preserve the Lorentz 
invariance, since the Hadamard regularization is performed in the Lorentzian rest frame of each of the 
particles [2.4] . However we have not been able to fix the kinetic ambiguity in the 3PN quadrupole moment 
using the extended Hadamard regularization. 
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in some subtle differences arising between the "case-by-case" computation of the elementary 
non-compact integrals in paper I and the systematic approach followed here which is based 
on the formula (|4.22jl . Of course, there is only one thing which is finally important, namely 
that these differences are completely encoded into some mere shifts of the values of the 
ambiguity parameters, see Eqs. (|5.7|) below.] 

The second source of ambiguity is "static" . It comes from the a priori unknown relation 
between the Hadamard regularization length scales, si and S2, and the ones, r[ and r' 2 , 
parametrizing the 3PN equations of motion in harmonic coordinates 0, |22| . The constants 
r[ and r' 2 come from the regularization of Poisson-type integrals in the computation of the 
equations of motion, and can be interpreted as some infinitesimal radial distances used as 
cut-offs when the field point tends to the singularities. Since we need the equations of 
motion when computing the time derivatives of the 3PN quadrupole moment, for instance 
in order to obtain the gravitational-wave flux, we must definitely know the relation between 
Si, S2 and r[, r' 2 . This relation constitutes a true physical undeterminacy within the various 
variants of Hadamard's regularization (either ordinary, pHS, hybrid or extended). Let us 
rewrite the R.H.S. of Eq. (J5.1j) by "artificially" introducing r[ and r 2 into the two slots 
of the pHS result. For doing this we use the known dependence of the pHS quadrupole in 
terms of the constants s\ and S2- This dependence is the same as in the case of the "hybrid" 
quadrupole (and indeed of any other of its regularization variants), and is given by Eq. 
(10.4) of paper I. Hence we have 

if^U ^ r } = ^^ In ( "J?) yfaf + 1 - 2 + • • • , (5.2) 



where a\ denotes the Newtonian acceleration and the dots indicate the terms that are inde- 
pendent of Si and s 2 (but which can depend on r ). Using this structure it is evident that 
the effect of changing si,«2 — ► r [, r 2 i n the pHS quadrupole is 

lff iS [s u s 2 ,r } = /^ HS K,r2,r ] 

+ 44 G^mf ^ fr\ \ ^ + x ^ ^ 



3 c 6 \sx 

We now argue, exactly like in Section X of paper I, that the most general admissible structure 
for the unknown logarithmic ratio in Eq. ()5.3|) is 

ln( r -±)=£ + K mi + m2 andl-2, (5.4) 
\sij mt 

where £ and k denote two new arbitrary ambiguity parameters, which are also a priori 
different from £ and k in paper I. The argument leading to Eq. ()5.4|) is essentially that 
the quadrupole moment should be a polynomial in the two masses mi and m 2 separately. 
Therefore, 

I,f S lH,S2,ro] = ifKSvro] 

3 c° \ mi J 

so we write the Hadamard-regularized 3PN quadrupole, depending on the three ambiguity 
parameters £, k and (, in the form 

Iij[si,a 2 ,r ;C] = I^ S [r'i,r' 2 ,r ] 
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+ 



44 G 2 m\ 
Y c 6 



m± 



l<->2. (5.6) 



We recall from paper I that, by contrast with the latter ambiguity parameters, the three 
scales r[, r' 2 and r are not physical and must disappear from the final results (when they 
are expressed in a coordinate-invariant way). 

Next let us compare the result, for general non-circular orbits, with the one of paper I 
which was obtained by means of the hybrid Hadamard-type regularization. 13 If everything 
is consistent, Eq. (j5.6J) should be in perfect agreement with paper I modulo a change of 
definition of the three ambiguity parameters, due to the use of the pHS regularization here 
instead of the hybrid regularization in paper I. We find that indeed there is a complete 
match for all the terms with those of paper I if and only if the ambiguity parameters £, k 
and ( are related to the corresponding ones £, k and ( in paper I by 

e = e + ^, (5.7a) 
k = k, (5.7b) 
C = C + ^. (5.7c) 

In view of the many differences between the present calculation and the one of paper I 
(e.g. in the definition of the regularization, the choice of elementary potentials, the way one 
computes non-compact support integrals), this agreement constitutes an important check 
of the lengthy algebra and the correctness of the result. In the following we prefer to come 
back to the original ambiguity parameters £, k and ( adopted in paper I, so we write the 
quadrupole moment as 



14 



+ 



V 

44 G 2 m{ f 1 mi+m 2 
T c 6 

l<->2. (5.8) 



Finally we present the result of the computation of all the terms in Eqs. (I3.4j) - (|3.6|) 
for 1 = 2 and general binary orbits. Unfortunately we find that the end expression of the 
quadrupole is very long in a general frame (with arbitrary origin), so we decide to present 
only the much shorter expression valid in the frame of the center of mass. The center-of- 
mass frame is defined by the nullity of the 3PN conserved integral of the center-of-mass 
vector deduced from the 3PN equations of motion in harmonic coordinates [3^|. For this 
calculation we use the relations between the general and center-of-mass frames given at 3PN 



For the purpose of the comparison we have redone the calculation of paper I in the case of non-circular 
orbits. In fact, although the end result of paper I is presented for circular orbits, most of the intermediate 
expressions in this paper are valid for general binary orbits. 

We employ the slightly abusive notation that Iij [s±, S2, ?"o; C] = hj r' 2 , tq\ £, k, C] when Eqs. (|5.4|l and 
(|57}l hold. 
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order in Ref. |47|. The structure of the 3PN center-of-mass quadrupole moment is 



A#i,r2,r Q ;£,/«,C] 



v m 



A 



24 v G 2 m 2 
7 c 5 r 2 



+2 



r r 24 v G m 



X(iXj) + B 

2 



XuV 



(5.9) 



Here we have explicitly displayed the "odd" 2.5PN radiation reaction contributions. The 
content of the "even" terms is given by the coefficients A, B and C, which generalize to 
non-circular orbits those given in Eqs. (1 1.3)— (1 1 .4) of paper I, and read 



A = 1 + 



+- 



29 
42 



29 v 
IT 



Gm 



Gm /2021 5947 



r 

Gm 

r 



V 756 

2 



V — 



756 
355 953 



5 
7 
4883 
756" 



+r 2 



253 
504 
Gm 



252 
1835 



126 



337 
252 



3545 2 
504 504 
131 907 1273 



1 

c 6 



r 

4561 



1 v 

756 756 



756 



7993 



11088 1584 



v + 



117067 



328663 



5544 



v 



11088 



v 



+v q 



+ 



Gm /307 94475 



; \ 77 

3 



Gm 



4158 
6285233 

+ 



v 



218411 
8316 



299857 



8316 



v 



207900 



44 

-_„(£ + 2k) 



+r 2 



Gm 



Gm 



2 

99 

2 



428 
105 
8539 
"20790 
1745 



34091 


3632 


V 




1386 


693 




44 


"(3 


i/ 




3 



13289 
8316 

r 



52153 
+ -^-v - 



2772 
187183 
83160 



4158 
16319 
5544 

605419 



4652 „ 54121 
v l - 



231 
311 



5544 



v 



16632 



v + 



99 
434909 



37369 



16632 



v 



2772 



2 Gm . 2 
+ it r 



757 5545 



5544 8316 



98311 
16632 



153407 
8316 



v 



B = v 

21 7 



(5.10a) 



15 Our notation is m = mi + m2 and v = mim2/m 2 ; x l = — an< l u * = dx l /dt = v\ — v\; v 2 = v 2 and 
r = n.v, where n = x/r and r = |x|; the STF projection is indicated by brackets surrounding the indices. 
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(5.10c) 



The 3PN quadrupole moment depends on £, k and C, on the constant scale r introduced into 
the general formalism defined for extended PN sources in Eq. (|2.7|) . and on the "logarithmic 
bary center" r' of the two Hadamard self-field regular ization scales r[ and r' 2 , defined by 



m In r' = m\ In r[ + m2 In r' 2 . 



(5.11) 



Unlike tq which cancels out in the complete waveform, already at the level of the general 
"fluid" formalism, and r' which represents some gauge constant devoid of physical meaning 
(see and paper I), the ambiguity parameters £, k and ( represent some genuine physical 
unknowns, which have recently been computed by means of dimensional regularization in 
Ref. (3^]. We shall show in the next Section that it is possible to determine a particular 
combination of these parameters in the context of Hadamard's regularization. 
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B. The 3PN mass dipole moment 



Recall from Section III CI that the mass-type dipole moment or "ADM dipole moment" 
Mi, which varies linearly with time (M, = 0), is the sum of two terms, 

Mi = Ii + 5Ii, (5.12) 

where J, is defined by the same general expression (|2.18aj) as for non-conserved moments 
but in which we set I — 1, and where the suplementary piece 5Ii is given by Eq. (j2.22b|) . 

We first concentrate our attention on the first part Li which is thus given, up to 3PN 
order, by the explicit expressions (|3.4|) - (j3.6j) with I — 1. We follow the same steps as for 
the quadrupole moment investigated in Section IV Al Repeating the arguments presented 
in Section X of paper I, we notice first that in the case of the dipole moment there is no 
ambiguous terms of the "kinetic" type. Actually one can easily check on dimensional grounds 
that the existence of such a term in /j, which would be proportional to v\ (plus 1 «-> 2), is 
impossible. Thus, unlike in the quadrupole case as shown in (|5.1|) . U is directly given by the 
result of the pHS regularization, 

h[si,s 2 ] = lP ns [ Sl ,s 2 ]. (5.13) 

Here s\ and s 2 are the two regularization scales coming from Eq. (|4.4|) . but as it turns out 
that there is no dependence on the cut-off scale ro in the dipolar case. To define the static 
ambiguity we must now re-express the dipole moment in terms of the particular equation- 
of-motion-related scales r[ and r' 2 . For this we use the dependence of the dipole moment in 
terms of the scales Si, s 2 , 

s 2 ) = ? ^ a{ In (^) + 1 ~ 2 + • - - , (5.14) 



3 c 6 1 \ Sl 

where the dots represent the terms that are independent of si and s 2 . Notice the factor 
22/3 instead of 44/3 in the quadrupolar case (|5.2|) . This yields immediately 

The ratio r[/si is a priori unknown but we remember that it has already served for the 
definition of two of our ambiguity parameters: £ and k, see Eq. (|5.4j) . Now we shall use 
in our present calculation of the dipole moment the same relation between r[ and Si as 
was used for the quadrupole moment. This means that we consider that the constants si 
and S2 parametrizing the Hadamard partie finie (|4.4|) have been chosen once and for all at 
the beginning of both our calculations of the quadrupole and dipole moments, where they 
take some definite meaning related for instance to the shape of the regularizing volumes B\ 
and B 2 which are initially excised around the two singularities when applying Hadamard's 
definition in the form of Eq. (J4.4j) . Thus we assume that s% and s 2 represent some unknown 
but fixed constants — having the same values for the two calculations of the quadrupole 
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and dipole moments. 16 When substituting the expression In (r[/ s\) = £ + k + k m^jmi into 
Eq. (J5.15)) we observe that the last term, which is proportional to the mass ratio m^jviix^ 
cancels out after applying the symmetry exchange 1 <-> 2. So we find that the dipole moment 
depends in fact on one and only one combination of ambiguity parameters, namely £ + k, 
and we obtain 

+ y^^^ + ^j + (5.16) 

Then we come back to the original definitions of paper I by using Eqs. (|5.7|) . and this leads 
to the following expression of the 3PN dipole moment: 

IiKriit + K] = L? HS [r[,r>] 

"3 V 22 



+ — — ^ (£ + «; + — J 4 + 1^2. (5.17) 



At this stage we have to worry about the extra contribution Sli present in Eq. (|5.12j) . 
From its expression given by (|2.22b|) we see that obtaining Sli at the 3PN order requires 
both E a and E a & with the full 3PN accuracy. By contrast, recall that the calculation of /; 
necessitated E at the 3PN order, but E a and E afe with only the 2PN and 1PN precisions 
respectively. Thus it seems that Sli cannot be obtained solely with the formulas developed 
in Section 11111 Notice that the expression of Sli involves an explicit factor B, and thus 
depends only on the presence of IR poles oc 1/B in the integrals. Consequently Sli can be 
computed by the same techniques as in Section IIVD1 i. e. in the form of surface integrals 
at infinity similar to Eqs. ()4.28|) or ([4. 30)1 . We have been able to prove that all the terms 
in Sli are separately zero up to the 3PN order. For all the terms we did know from using 
the results of Section ITTT1 we have made a complete calculation, and for the other terms we 
looked at their allowed structure in terms of the basic potentials V, Vi, Wij, • ■ ■ , invoking 
dimensionality arguments but leaving aside the unimportant numerical coefficients, which 
was sufficient to check that the corresponding surface integrals are exactly zero for all the 
terms. Thus, we conclude that Sli = at 3PN order, hence 

M i = I l + 0(c- 7 ) , (5.18) 

which finally results, from the detailed evaluation of all the terms in Eqs. (I3.4|) - (j3.6|) . in 

Mi = m iy \ 



6 We tried to further extend this type of argument to the constants s\ and S2 which were used in the 3PN 
equations of motion 0, ■ This implied that we had to use for the wave generation the same regular- 
ization as in the equations of motion, i.e. the extended Hadamard regularization |2.'{ll'2 1( . Unfortunately 
this program, whose aim would have been to determine all the ambiguity parameters within Hadamard's 
regularization (£, n, ( and also A), did not fully succeed. Nevertheless a less ambitious part of the program 
did succeed, and this is what we show here. 

7 The two masses mi and m-i are located at the positions yi and y2, the unit vector between them is 
n i2 = (yi — y2)/^i2 with ri2 = |yi — ya], the two coordinate velocities are vi = dy±/dt and V2 = dy^ldt^ 
and V12 = Vi — V2. Euclidean scalar products are denoted by parenthesis, e.g. (jiyivx) = ni2.vi. 
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The case of the conserved dipole moment is interesting because it offers us a very good 
check of the calculations. Indeed let us argue that Mi, which represents the distribution 
of positions of particles as weighted by their gravitational masses m g , must be identical to 
the position of the center of mass Gi of the system of particles (per unit of total mass), 
because the center of mass Gi represents in fact the same quantity as the dipole Jj but 
corresponding to the inertial masses ni{ of the particles. The equality between mass dipole 
Mi and center-of-mass position Gi can thus be seen as a consequence of the equivalence 
principle = m g , which is surely incorporated in our model of point particles. Now the 
center of mass Gi is already known at the 3PN order for point particle binaries, as one of the 
conserved integrals of the 3PN motion in harmonic coordinates (we neglect the radiation- 
reaction term at 2.5PN order). We recall that Gi, given explicitly in Ref. [35j], depends 
on the regularization length scales r[ and r' 2 which are the same as in the 3PN equations 
of motion 21, 2^1 and therefore the same as in our result (]5.19 | ) - - by definition of the 



ambiguity parameters £ and k. However Gi was found in Ref. [35| to be free of ambiguities; 
for instance the ambiguity parameter A in the 3PN equations of motion disappears from the 
expression of Gi. Let us therefore impose the equivalence between Mj and G i} which means 
we make the complete identification 

MMSvZ + k] =Gi[riy 2 ], (5.20) 

in which we insist that the constants r[ and r' 2 appearing in both sides of this equation are 
the same. Comparing Mj with the expression of Gi given by Eq. (4.5) in |35|, we find that 
these constants r[ and r' 2 cancel out, and that Eq. ()5.20|) is verified for all the terms if and 
only if the particular combination of ambiguity parameters £ + k is fixed to the unique value 

(5-21) 

^ 9240 V ; 

This result is obtained within Hadamard's regularization. It shows that, although as we have 
seen Hadamard's regularization is "physically incomplete" (at 3PN order), it can nevertheless 
be partially completed by invoking some external physical arguments — in the present case 
the equivalence between mass dipole and center-of-mass position. 

More importantly, we find that Eq. 1)5.21)1 is nicely consistent with the calculation of the 
ambiguity parameters by means of dimensional regularization whose results have 

been given in Eq. (jl.lj) . The dimensional regularization is complete; it does not need to 
invoke any "external" physical argument in order to determine the value of all the ambiguity 
parameters. Nevertheless, it remains that our result (|5.21|) . based simply on a consistency 
argument (within the over-all scheme) between the 3PN equations of motion on the one hand 
and the 3PN radiation field on the other hand, does provide a verification of the consistency 
of dimensional regularization itself. 
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